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The quantum algorithm of Kieu 
does not solve the Hilbert's tenth problem 



Boris Tsirelson 



Abstract 



Recently T. Kieu [jjl claimed a quantum algorithm computing some 
^ I functions beyond the Church- Turing class. He notes that "it is in fact 

' widely believed that quantum computation cannot offer anything new 

^ . about computability" and claims the opposite. However, his quantum 

algorithm does not work, which is the point of my short note. I still 
believe that quantum computation leads to new complexity but retains 
the old computability. 



1 The algorithm of Kieu 



^ ' Obtaining the ground state of a Hamiltonian could be treated as a basic 

D^. operation (feasible by definition), but such an algorithm should be called 

algorithm with a ground-state oracle rather than quantum algorithm. By the 
'O I quantum algorithm of Kieu I mean the method described in Sect. "Adiabatic 

; evolution" of §.[1 

Here is the relevant portion of the algorithm. In order to find the ground 
state 1^') of a given Hamiltonian Hp, one uses the adiabatic evolution from 
the known ground state 1^'/) of an initial Hamiltonian Hj: 



\g) = ^lim r exp\^-i Hit) dt j \gj) 

up to a phase; that is, a state vector \g{t)) evolves by the Schrodinger 
equation 

'git)) =-tn{t)\g{t)) forO<t<T, 



dt 



W)) = \9i) , 



^Though, he did not specify the initial Hamiltonian Hi. We'll see that his algorithm 
fails for every Hj. 
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where 

Hit) 

giving \g{T)) ^ \g) 

2 The goal of the algorithm 

The quantum algorithm of Kieu is intended for finding the global minimum 
of a function^ P{x,y,z) of, say, three variables x,y,z each running over 
nonnegative integers {0, 1, 2, . . . }.0 To this end, the adiabatic evolution finds 
the ground state of the Hamiltonian Hp whose eigenvalues are the numbers 
P{x,y,z). 

Let us use a single- variable function P{x), just for simplifying notation. 
Our Hilbert space is spanned by basis vectors |a;) for x G {0, 1,2,...}, and 
Hp is a diagonal operator. 

Hp \x) = P{x) \x) for X e {0, 1, 2, . . . } . 

The user of the algorithm has no upper bound for the minimizer Xmin 
of the function P. You see, if an upper bound is known, then the global 
minimum can be found classically, just by computing and comparing a finite 
number of values P{x). The goal of the algorithm is finding the global 
minimum 

Pixmin) = min P{x) 
xe{o,i,2,...} 

and the minimizer x^i^, during a time T that does not depend on P. That is 
crucial: the whole infinite domain {0, 1,2,...} must be effectively examined 
during the finite time T. According to Kieu, that is possible due to infinite 
dimension of the Hilbert space. 

If Hi is also diagonal in the given basis ( \x) )x=o,i,..., then clearly \g) = 
\gi) (up to a phase); nothing happens. This is why operators Hj and Hp 
must be non-commutative, which is noted by Kieu, Sect. "Gap estimation" 
of [||]. Matrix elements of Hj must connect basis vectors. However, they 
cannot effectively connect first basis vectors with remote basis vectors \x) 
for arbitrarily large x, as one could guess. That is the obstacle. See the next 
section for a proof. 

^Kieu considers only polynomials P, which is essential when implementing the Hamil- 
tonian Hp, but irrelevant when finding the minimum by adiabatic evolution. 

^For example, the function P{x,y, z) — [{x + l)-^ + (y + 1)'^ ~ (-^ + 1)'^)^ is evidently 
related to the Fermat's last theorem. 



l-^Hj + ^Hp forO<t<T, 

7^(0) = Hi , n{T) = Hp , 
(up to a phase) provided that T is large enough. 
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3 The goal is not achieved 

Consider such a function -P(-):0 

P{x) 



1 if X 3^miru 



otherwise. 

In order to disclose its minimum at Xmin (unknown to the user), the adiabatic 
evolution is harnessed: 

^\g{t)) =-tn{t)\g{t)) forO<t<T, 



dt 



W)) = \9i) , 
nit)=(l-^)Hr + ^Hp. 



Compare \g{-)) with the solution \go{-)) of a simpler equation 
di 



Note that \go{t)) = \gi) up to a phase (varying in time), since \gi) is an 
eigenvector for Hj, therefore for T-Co{t). We have 

^ \9{t) - go{t)) = -tUit) \g{t) - go{t)) + i^Hp \go{t)) ; 

^11 \g{t) - go{t)) f = 2 Re ( g'{t) - g',{t) \ g{t) - go{t) ) = 

= -2 Im ( ^Hpgoit) \ git) - goit) ) <2\\Hp \goit)) \\ ■ \\ \g{t) - g^{t)) \\ ; 

^11 |^(t) - g^{t)) II < \\Ep \g^{t)) II = \\Ep \gj) \\ ; 

\\\g{t)) - \g^{t))\\ <T-\\Hpgi\\. 

The arbitrary parameter Xmin, unknown to the user, does not influence T, 
Hi and gj. However, 

\9l) II = I ( Xmin \gi) \ ^ . 

^'min 

For large Xmin we see that \g{T)) is close to \gQ{T)) ; the latter, being equal 
to \gi) up to a phase, cannot be close to |xjnin) • The adiabatic evolution 
fails to disclose the minimum. 



■*Kieu considers only nonnegative functions. Well, my argument works also for the 
nonnegative function P{x) + 1. 
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